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We study the heavy-light mesons in a relativistic model, which is derived from the Bethe-Salpeter 
equation by applying the Foldy-Wouthuysen transformation on the heavy quark. The kernel we 
choose is based on scalar confining and vector Coulomb potentials, the transverse interaction of the 
gluon exchange is also taken into account in this model. The spectra and wave functions of D, D s , 
B, B s meson states are obtained. The spectra are calculated up to the order of 1 /mQ, and wave 
functions are treated to leading order. 
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I Introduction 


Recently, great progress has been made in the mea¬ 
surement of the spectrum of heavy-light quark-antiquark 
system, especially in the charm sector. Several new ex¬ 
cited states were established in both D and D s mesons 
M- Many theoretical models were propounded to iden¬ 
ti fy th e new discovered states and explain their properties 
[9H14j. The spectroscopy provides a powerful test of the 
theoretical predictions based on the quark model in the 
standard model. 

Heavy-light quark-antiquark system Qq plays an im¬ 
portant role in understanding the strong interactions be¬ 
tween quark and antiquark. Double-heavy meson sec¬ 
tor is rather mature, and nonrelativistic potential models 
have been proven extremely successful for the description 
for heavy quarkonia 0. As for heavy-light meson sys¬ 
tem, one needs a model that can include the relativistic 
effects of the light quark. 

In order to describe the relativistic effects in the bound 
states, we resort to the Bethe-Salpeter equation. Many 
studies of mesons were carried out in the Bethe-Salpeter 
approach 0-0- It is very difficult to solve the Bethe- 
Salpeter equation, especially when we consider states 
with large total angular momentum number J. We can 
reduce the Bethe-Salpeter equation to an approximate 
but easier form, so as to study mesons systematically. 
We apply the Foldy-Wouthuysen transformation on the 
heavy quark and expand the interaction terms to order 
1 /toq. The instantaneous approximation can only con¬ 
tribute to corrections of order (1 /toq) 2 if one considers 
the heavy quark on shell as in the heavy quark effective 
theory (HQET). In this paper we take the instantaneous 
approximation for the Bethe-Salpeter equation. 

One can bring the instantaneous Bethe-Salpeter equa¬ 
tion into the Breit form by taking further approxima¬ 
tion. The Breit interaction is widely used to study quark- 
antiquark systems as an nonrelativistic or semirelativistic 
model ['2 1 h 23I . While in this paper, we take the instanta¬ 
neous Bethe-Salpeter equation in its original form, from 
which one can derive a more reliable model that preserves 
the relativistic effects of the light quark. 


We treat the equation in the 1 /toq expansion to the 
leading order. The parameters in the equations are de¬ 
termined by fitting the spectra of the heavy-light meson 
states presented in Particle Data Group (PDG) [24j]. The 
identifications of the newly observed heavy-light meson 
states are given based on the theoretical prediction. 

The organization of this paper is as follows. Section II 
is for the relativistic model and the effective Hamiltonian 
of the heavy-light quark-antiquark system. In section III, 
we solve the relativistic wave equation. Section IV is for 
the numerical result and discussion. In section V we have 
a brief summary. 


II The model 


The Bethe-Salpeter Equation for the quark-antiquark 
system can be written as 25J: 


■ dV 
(2ny 


(fti - nn)x(p)(h + m 2 ) = 


-K{p,p',P)x{p'), (1) 


where the total momentum P (in the center-of-mass sys¬ 
tem, P = (E , 0 )) and the relative momentum p are de¬ 
fined as: 


P = Pi +P2, 

P = aip 2 - a 2 pi 


with 


TO,; 


toi + m-2' 


,* = 1 , 2 . 


( 2 ) 

( 3 ) 

( 4 ) 


It is necessary to assume some form for the kernel in 
order to reduce the Bethe-Salpeter equation. Here we 
take the kernel as the simplest form [26J, which can be 
reduced to that used in Ref. 0. The kernel can be 
written as: 


K(p,p',P) = 7 (1) • 7 {2) V v (k 2 ) + V s (k 2 ), (5) 


where k is the transferred four-momentum, 
k = p — p'. 


( 6 ) 
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Inserting Eq. (0 into Eq. 0, the Bethe-Salpeter equa¬ 
tion can be written as: 

(j/>i - wi)x(p )(^2 + m 2 ) = 


(27T) 4 


V v {k 2 ) 7 (1) • x(p + /c)7 (2) + V s (k 2 )x(p + k) 


(7) 


The above equation can be transformed to a more con¬ 
venient form that is symmetric with respect to particles 
and antiparticles by using the charge conjugation trans¬ 
formation. We rewrite Eq. 0 as: 

{i>\ - mi)(/S 2 + m 2 ) T x(p) = 


' d A k 
(2^ 


V v (k 2 b (1) • 7 {2)T \{p + k ) + V s (k 2 )x(p + k ) 


After applying the condition 

C^ T C~ l = - 7 ", 

we arrive at 

(A - mi)(^2 - m 2 )i/j(p) = 
f d A k 


( 8 ) 


(9) 


(27T) 4 


V v (k 2 ) 7W • 7 ^ip{p + k) — Vs(k 2 )i/j(p + fc) , 

( 10 ) 


where ^(p) is the charge conjugated wave function, and 
the minus sign appeared before 14 indicates that the con¬ 
finement force is always attractive. 

By taking the instantaneous approximation, we can 
perform the p° integration and decompose the instanta¬ 
neous Bethe-Salpeter equation into four equations. De¬ 
tails of this derivation can be found in Refs. dum. 

The projection operators are defined as 


A« = -(1 ±h i ), 


hi = 


H, 


OJi 


Ui = \ p 2 


i = 1,2 


with 


H\(p) = —a^-p + p^mi, 
H 2 {p) = a*- 2 ' 1 • p + /3^m 2 . 


( 11 ) 

( 12 ) 

(13) 


(14) 

(15) 


Applying the projection operators, the four coupled equa¬ 
tion can be written as 


(E-wi- w 2 )<j> ++ (p) 

= -2mY ++ (p), 

(16) 

(E + uji + u 2 )<j> — (p) 

= 27riT__(p), 

(17) 

<P+-(P) 

= o, 

(18) 

<t>-+{p) 

= o, 

(19) 


where 




dp°ip(p°,p), 


( 20 ) 


T(P) = 


(2?r)4 7 ° 7o 


K(-fe 2 ) 7 (1) -7 


( 2 ) 


-K(-fe 2 )] <Kp + fc), 

<f>±± = A^A^V, r ±± = A«A£ 2) r. 


( 21 ) 

( 22 ) 


One can combine the four coupled equations and ob¬ 
tain the equation of (f> instead of its projective compo¬ 
nents <j)±± as in Ref. (30], that is 


(E - Hi(p) - H 2 (p))<f>[p) = — 27r*A(p)r(p), (23) 


where 


A (p) = A«Af-A w A< 2) 

= \{hi{p) + h 2 {p)). (24) 

If one makes the approximation 

A(p) -)■ 1, (25) 

Eq. (l23l) can be bought into the Breit form. The Breit 
equation can be used as an nonrelativistic or semirela- 
tivistic approach to study the interactions of quarks. In 
order to take the relativistic corrections of the light quark 
into account, we keep A (p) in its original form instead of 
taking it as 1. 

Eq. (l23l) can be transformed into coordinate space 
according to 

4>(r) = j e^(p)- (26) 

The coordinate form of Eq. (E^l) is 

(Hi + H 2 + ^(h 1 + h 2 )U^j 4>{r) = E<f>(r) (27) 

with 

U{r) = U\{r) + U 2 (r), (28) 

Ui(r) = K(r)+/? (1) /3 (2 V s (r), (29) 

U 2 (r) = — + (cc (1) • r)(a (2) ■ r)]V v (r), 

(30) 

where V v ( r ) and 14 (v) are the Fourier transformations of 
V v (-k 2 ) and I4(—fe 2 ), respectively. 

The wave function <j)(r) has two spinor indices, i.e. 
4x4 = 16 components, it is difficult to solve an 
eigenequation with so many components. For the heavy- 
light system, one can reduce the eigenequation by apply¬ 
ing the nonrelativistic approximation of the heavy quark 
Q and expand the interactions in the order of 1/mg. 
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In an elegant and systematic way the reduction can be 
achieved by using the Foldy-Wouthuysen transformation. 
Here we designate the heavy quark Q and the light an¬ 
tiquark q by the index 1 and 2 in Eq. 03, respectively. 
Thus the transformation is performed on the heavy quark 
with the superscript “1”. We reduce the Hamiltonian in 
the Dirac representation, where 

“=(„ ")'' s =( 1 = *)■ (31) 

If the original Hamiltonian is written in the form 


If the above equation is true, then 

(hi - h 2 )<\> = i(/i! - h 2 )(h 1 + h 2 )p = 0, (39) 

on the other hand, if Eq. ((371) is true, we take ip = hi<p, 
then 

-(hi + h 2 )p = i(l + hih 2 )(t) = 4>. (40) 

Applying the equivalent form of Ea. (l37l) . Eas. (l36l) can 
be transformed to a new equivalent equation 


H = Pm + £ + 0, (32) 

where O is the “odd” operator, typical examples are the 
matrices a and 7 , while £ is the “even” operator, exam¬ 
ples of this class of operators are 1, /3 and S. Accord¬ 
ing to the Foldy-Wouthuysen transformation, the trans¬ 
formed Hamiltonian reads 


H 


u- [ hu f 

^ + £ + A 03 + _L| [cw] , 0] 


8 m 3 


O 4 


(33) 


The Hamiltonian in Eq. m is 

H = Hi + H 2 + ^(hi + h 2 )U. (34) 

If we perform the Foldy-Wouthuysen transformation di¬ 
rectly on Eq. (1341) . the result can be awkward. What we 
want to achieve is to obtain a reduced and simple form 
of the Hamiltonian for the heavy-light quark-antiquark 
system. Here we show that it is possible. It should be no¬ 
ticed that Eq. (l27l) is not equivalent to the instantaneous 
Betlie-Salpeter equation, i.e. the four coupled equations 
Eqs. ©-03), as pointed out in Refs. 3.1. [32|. We 
find that our goal can be achieved by employing the con¬ 
straint of Eqs. m and (fl9l) which are dropped in Eq. 

m- 

Subtracting Eq. (1191) from Eq. urn we get 

(hi — h 2 )(j) = 0. (35) 

We write Eq. (l27l) together with Eq. (l35l) as: 


(^Hi +H 2 + i(hi + h 2 )U - E^J i(hi + h 2 )ip = 0, (41) 

it is easy to verify that the above equation is equivalent 
to the two coupled equations (l36l) and (l37l) . Furthermore, 
it is equivalent to 

^—(1 + hih 2 )(uji + uj 2 ) + — (hi + h 2 )U — (hi + h 2 ) 
~(hi + h 2 )E S jip = 0, (42) 

Noticing the relations 


(1 + hih 2 ) — (hi + h 2 ) 

— hi + ii2i 

(43) 

(1 + hih 2 )hi 

= hi + / 12 , 

(44) 

(1 + hih 2 )h 2 

= hi -j- / 12 , 

(45) 


we can have two different equivalent forms of Eq. (1421) , 
they are 

— (1 + hih 2 ) + ui 2 + — (hi — (hi + h 2 ) 

"(hi + h 2 )E^j ip = 0, (46) 

and 

+uj 2 + i(/ii + h 2 )U^(hi + h 2 ) 

"(hi + h 2 )E\ i(l + hih 2 )ip = 0. (47) 


^H 1 +H 2 + ^(hi + h 2 )U - E^j <j>(r) 
(hi - h 2 )4>(r) 


0, (36) 

0. (37) 


Here we guess a form that is equivalent to Eq. (l46l) 
and Eq. (1471) by employing the common part of them, 
that is 


It is easy to verify that the above two equations are 
equivalent to the instantaneous Bethe-Salpeter equation, 
i.e. the four coupled equations Eqs. Here we 

expect to obtain one single equation that is equivalent to 
the above two equations by combining them. 

Eq. (1371) is equivalent to 

4> = 7^(hi + h 2 )ip. (38) 


+ oj 2 + —(hi + h 2 )U -(hi + h 2 ) 

~\(hi + h 2 )E^j -tp = 0. (48) 

Left multiplying the above equation by ^(1 + hih 2 ), 
we get the form of Eq. (147)1) . On the other hand, we can 
obtain the form of the above equation from Eq. fl?D if 
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we take ip = ^(1 + hih 2 )ip. That is to say, Eq. (l48l) is 
equivalent to Eq. (1461) and Eq. (l47l) . thus equivalent to 
Eqs. dH and @3, i.e. the instantaneous Bethe-Salpeter 
equation. 

By using Eq.fHJ and Ea. diSll . we can obtain two other 
equivalent equations, respectively. 

bq + u>2 + — (Iii ^ 2 )^ 2 (^H "T h 2 ) — hiE^j ip = 0, 

(49) 

+w 2 + 7;{hi + h 2 )U-(hi + h 2 ) - h 2 E^j ip = 0. 

(50) 

In the foregoing paragraphs, equivalent equations of 
the instantaneous Bethe-Salpeter equation are obtained. 
In Appendix A, the reduction of the instantaneous Bethe- 
Salpeter equation for nonrelativistic systems is discussed. 
We find that the reductions of the Hamiltonian with the 
approximation (1251) . i.e. the Breit interaction, and with¬ 
out the approximation (1251) . i.e. the choice in this paper, 
do not differ until order 1 /miirij. While for the case of 
heavy-light system we focus in this paper, the reduction 
results of the two scheme can be essentially different. 

Now we return to the reduction of the Hamiltonian for 
the heavy-light system. Eq. (l49l) can be written as (the 
wave function in the equation is omitted): 

— hiE +a w i + w 2 + -(^i + h 2 )U — {hi + h 2 ) = —in 1 , 

(51) 

where 

awi = oji — mi, (52) 

now we consider performing the Foldy-Wouthuysen 
transformation on Eq. m, the left of it can be rewritten 
as: 

-f,m E -E - 1 ) pm - e A AAA 

+A£Ui + U} 2 + ^{hx + h 2 )U^{hi + h 2 ) 

= -0O-)E + £ + O, (53) 

the odd and even operators in the above equation are: 

£ = -E - 1^ +AW1 + 0J 2 



and 



where h.c. stands for Hermitian conjugate. 

By using the new form of Eq. m 

H = -/?«£ + £ - ^O 2 + ^[[0,£\,0\ + ■■■ (56) 

we can expand Eq. m to order 1/E 

-p {1) E + Hq + Hx +# 2 = 0, (57) 


where 




Ho = 

Wl + UJ 2 + 1 

0«+/i 2 )c7i^(/ 

3 (1) + h, 

Hi = 

1 / 

I U 2 \ (/3 (1) + h 2 ) 

I -1- h.c. 


2 \ 0 J 1 ) 

1 2 V ) 



I—{« (1) ' h/? (1) + h 2 )U 2 ~(/? (1) + h 2 )}, 

Ami A A 


H, = - 


0W 

2 E 


i(/3W+/i 2 )t7 2 l(/3 (1) +/i2 


(59) 

(60) 


The last term H 2 is of the order Vf. Comparing to the 
other terms, its correction is suppressed by the coupling 
constant, here we omit it. Considering the substitution 
pw —> 1 in the Foldy-Wouthuysen transformation, Hi 
can be reduced to 

If a^ 1 ) • p\ 1 

Hi =► — ^ — ju 2 -{l + h 2 ) + h.c. 

1 1 / Ot^ • p\ 1 

+ - 2 (-i-hh 2 )-(-^)u 2 -{i + h2 ) + h,. 

1 1 / QjW • T)\ 1 

= 2< I+ «2(-^r) f/2 2 (1+ « + '“- 

= |(l + 'l2)k-^- £ ,C2(hl + /l2), (61) 

Z Z 777-1 z 

then H is transformed to 

H = CJi + 072 + — (1 + ^ 2 ) f/i— (1 + ^ 2 ) 

1 1 • V 1 

+ ^(1 + h 2 )-{ - -,U 2 } ~{l + h 2 ), (62) 

Z Z 777-1 Z 
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and we obtain the final Hamiltonian 


H = H 0 + H' (63) 

H 0 = uJi+uj 2 + ^(l + h 2 )(V v +/3^ 2) V s )^(l + h 2 ) (64) 


H ' - 2 ( 1 + /l2 ){ 4mi [(a (2) + (« (2) 

■ r)r)■p 

+ icr 1 ' 1 ' 1 ■ (ot { - 2 ' > + (a*- 2 ^ • r)r) x p 

+ h.c.} — (l + h 2 ) 


(65) 

The scalar and vector potentials we choose have the 

simple form 


V a (r) = br + c, 

(66) 

t/ i \ A(Xs 

VJr) = -—' 

(67) 


The potentials are chosen to have a Coulombic behavior 
at short distance and a linear confining behavior at long 
distance. 


Ill Solution of the Wave Equation 


In this section, we solve the eigenequation of Hq and 
then deal with the perturbative corrections of H'. 

In order to generalize the algorithm which was used in 
Refs. [H, 0, we prove a relation in the first place. 

Let f l(p) be the pseudo-differential operator function 
and f l(k) be the normal function, here p and k stand for 
the modules of momentum operator p and momentum 
k, respectively. Let 4’( r ) be a function of r that can be 
written as 


ip(r) = 4>(r)Yi m (r), 

with the help of 

A 3 !' r — r '\ — f ^ ^ ik-(r-r') 

1 V ~ J (2ir) 3 3 


we have 


( 68 ) 


(69) 


f l(p)ijj(r) = Q(p) J d 3 r' S 3 (r — r')ijj(r') 


= J dV / ^3 (70) 

The exponential factor e lk r can be decomposed into 
series of spherical harmonics 

e lkr = 4t TY l i l 3 i(kr)Yt n (k)Y lm (r), (71) 


where ji is the l -th order spherical Bessel function, Y) m (r) 
is the spherical harmonics, which satisfies the normaliza¬ 
tion condition 


(f') — di 1 i 2 S mirn2 , 


(72) 


with r the unit vector along the direction of r. 

Inserting Eqs. (168)) and HID into Eq. m, after sim¬ 
plification we arrive at 

Vt{p)<t>(r)Yi m (r) 

= ~ J dr ' r ' 2 J dk k2 tt( k )ji( kr )ji{ k r')<i){r')Y lm (r). 

(73) 

If we take fl(p) = 1, then O(fc) = 1, the above equation 
is simplified as: 

cj)(r) = — J dr'r' 2 J dkk 2 ji(kr)ji(kr')(j)(r'). (74) 

The above equation is true for arbitrary function <fi(r ), 
which gives the orthogonality condition for the spherical 
Bessel function 

~ r ' 2 J dkk 2 ji(kr)ji(kr') = 6(r — r'). (75) 

Furthermore, if we take </>(r) = ji(k'r), we can transform 
Eq. ( 1731 ) to a simple form by using Eq. 

ft(p)jl(k'r)Yi m (r) 

= j dkil{k)jt(kr)Y lTn (r)6{k - k'). 

= n(k')ji(k'r)Y lm (r) (76) 


that is 


tt(p)ji(kr)Yi m (r) = n{k)ji{kr)Yi m {r). (77) 


The above equation is actually an eigenequation of fl(p). 
If we expand the unsolved wave function in terms of the 
spherical Bessel function, the above relation can do a 
great help in solving the wave equation. 

Now we turn to solving the eigenequation of Hq. The 
spinor part of Eq. (11071) can be rewritten as 


4/(r) 


Ik A 


r )y7l 4 ^ 
*/(»■) ) ’ 


(78) 


by using the formula 


er • p = i(cr ■ r) f ^ (79) 

\ r dr J 

we have 

(a- ■ p)g(r) y™l A = i g(r ) y™* B , (80) 

(<7 • p)f(r) y™J B = -i g{r) y™’ A ,{ 81) 


where the relations 


(<* 


~ \ rrij rr 

r ) Vj,L = 


3,1b 


m j 

r > Vj,i B : -y 


3,1a 


(82) 

(83) 
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are used. 

Inserting Eq. da and m into Eq. (1641) . we can 
rewrite Hq in the matrix form 

H o=u 1 +uj2 + -(l + (V v + /3 (2 V s ) (l + , 

(84) 


with a n the n-th root of the spherical Bessel function 
ji{x ) = 0. N n is added in accordance with the normal¬ 
ization condition in Eq. (11111) . In the limited space, the 
momentum k is discrete, we can have the relevance 



where 



m2 I i & P 

U12 U>2 

&P m 2 , ^ 

Ul 2 U>2 ' 


V v + p^V s = 




then Hq can be written as: 


ffo=(“ 1+ “ 2 , ) 

V Wi + w 2 ) 

1 (H a H b \ 

4 \H C H d )' 


(85) 

( 86 ) 


(87) 


where 

H a = (— + 1)(K+K)(— +1) 

+ — (K-K) —, (88) 

UJ2 W2 


H b = (— + i)(k„ + v;)— 

W2 W2 


+ 

(T V, 

— ~(V v 

-Vs)(- 

-^ + 1), 

(89) 


UJ 2 

w 2 

H c = 

Hi 



(90) 

H d = 

(-— + 

i)(K - 

ViK-^ + D 



LU 2 



+ 

a p (V v 

+ V 8 )‘ T 

■p 

(91) 


W2 W2 


Since the quark and antiquark are bound in the me¬ 
son, when the distance between them is large enough, 
the wave function drops dramatically. The wave func¬ 
tion will effectively vanish at a typically large distance 
L , the quark and antiquark can be viewed as being re¬ 
stricted in a limited space, 0 < r < L. The unsolved 
functions f(r) and g(r) can be expanded in terms of the 
spherical Bessel function in the limited space, we write 
the basis as 




—Ji (^) 

N a M a[ l ’ 
i ■ 

N B^ B< y L ) 



VHb 


(92) 

(93) 


where the superscripts “A” and “ B” stand for the upper 
and lower parts of Eq. (j78j) , and N n is the module of the 
spherical Bessel function 


m = 


, / /2 ■ / u n' \ 2 

drr ji {——) , 


(94) 


The eigenfunction of Hq can be expanded in the or- 
thonormalized basis 

oo oo 

^ 9i^i + fai’a ■ (96) 

2=1 Oi=l 

In the numerical calculation the above summation can 
be truncated at a large integer N, 

N N 

2=1 06 = 1 

According to Eq. (HB, we can rewrite the eigenequa- 
tion of Hq in the representation of {ipf,ipa}- In this 
representation, the operator Hq has its matrix form 


Hq 


( < U> 1 + 0 l>2 >ij 

\ < U>i + W2 >a/3 

1 f ^ H a O ij < H b > i/3 \ 

4 \ < H c > a j < Hd O a /3 J 


(98) 


Applying Eqs. (1771) . (15U1) and (I5T1) with the normaliza¬ 
tion condition Eq. (II 101) , we can get the matrix elements 
of Hq easily. 


< Uli + W2 >ij — 

< U>\ + UJ2 >aB = 


<*(■£) + «*(■£) 


S ijt (99) 




1 


Wl( ^ )+W2( ^ } 


S a p, (100) 


here we define a symbolic notation 

rL 


(<t>(r))m, lA ;n,i B =l dr r*j lA (^f), 


( 101 ) 


then we have 


^ Hn Oil'- 


m 2 


a 'v- N a n a I (°£ 

1 o W(-^, 

X ^ V v + V s O ,1a 

1 1 1 

XI i — + i)' l - v -- v -> 


+ 1 


m 2 


(<■ 

uj 2 {- 


+ 1 


k + 1 d 

r dr 


i-ilA. 

( 102 ) 
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< H b >ip= 


1 


7712 


N±N§ 


+ 1 


1 


X ( (V v + Vs) 


p W a -t) ) utft) 

k — 1 d 
r dr, 


i, l a; PJb 


1 


1 


-7712 


J'W ,*,(!£) y nC _L } 


+i 


H c > aj — ^ H\) !> j ay 




(103) 

(104) 


In Eq. (11081) . y ; m±1 / 2 j s the spherical harmonics. For 
j = l ± 1/2, we have 

fdn(y%)iy% = 1 . ( 110 ) 


Thus the normalization condition of Eq. (|107|) is read 


as 


poo 

/ ^TT 2 +9n,l,j) = !■ (Ill) 

Jo 

For a state with quantum number j, the operator A' 
can have two opposite eigenvalues ±(j + 1/2) with 


< >«/3= ( -S- +1) ( —^ +1 


^ ^ Vy ^ocjIb 

1 1 1 

n b n b i a f \ TofT 

JV « iV /3 w 2 (-f-) W2 (-#-) 




" fc — 1 d 


r dr 


-31 (^ + V 8 ) 


fc — 1 d 


r dr, 


Diagonalizing the Hermitian matrix of Hq , we can get 
the eigenenergy of Hq and the coefficients gi, fon which 
are defined in Eq. ©a, then the eigenequation associated 
with H 0 is solved and the eigenfunction is obtained. 

Now we discuss the perturbative corrections of H'. It 
is easy to verify that the operators 

{j ) jz, K, S/} 

are a set of mutually commuting operators which com¬ 
mute with Ho, where j = L + S '* 2 ), S * 2 ^ = i£* 2 \ 
K = /3(2)(S (2) • L + 1), S = ler«. Then the eigen¬ 
states of Hq can be labeled by the corresponding set of 
quantum numbers {n,j, rrij, k, s} and the eigenequation 
associated with Hq can be written as 


M 


(o) 


o = E ni^Sk ( 106 ) 


where 


>*• (107) 


with 


( k + Y m ~ 1 / 2 \ 

3 ' 1 V k 7,l,mYr +1/2 J ’ 


U± — 


+V / ^^^' = Z + 1 /2 

3=l~ 1/2- 


-F 


7=Fm+l/2 
21 + 1 > 


(108) 


(109) 


l = 


j + 1/2, k = +(j + 1/2) 
j ~ 1/2, k = —(j + 1/2), 

where the quantum number k is defined as 


AN/ 


(°) (r) = —k'I> [ ■ u, (r) 

n,k,j,rrij / n,k,j,rrij ,s\ /* 


( 0 ) 


( 112 ) 


(113) 


cx,Ib',(3,Ib 

(105) 


The zeroth order invariant mass E® k • is equivalently 
determined by n,j,l , and the eigenstate has the parity 
P = (—1) I+1 . The Solution of the eigenequation associ¬ 
ated with Hq is detailed in the next section. 

The perturbative term H' does not commute with any 
of the operators introduced in solving the eigenequation 
of H 0 , but still commutes with 

{J 2 ,J Z ,V} 

where J = j + S and V is the parity operator. Thus the 
eigenstates of the total Hamiltonian H = Hq + H' can 
be labeled by the set of quantum numbers {n, J, Mj, P}. 

With the help of Clebsch-Gordan coefficients, we can 
compose the basis states with the quantum numbers set 
{n, k,j, J, Mj} by combining the eigenstates of Hq, 

" E <W/2(H4) 

rrij ,s 

then calculate the corrections and mixings in the ob¬ 
tained basis. 

The correction in first order perturbation can be writ¬ 
ten as 


? ( 0) 


E n , u ,j = Ef”, , + —SE. 


1 -y(l) 


m i 


(115) 


Mixing can happen between states with same quantum 
numbers J , Mj , P, but with different quantum number j. 
In this paper the mixing between states with j = l± 1/2 is 
considered. Mixing can be described by the mass matrix. 
The mass matrix is calculated perturbatively in the basis, 
which can be written as 


H = 


H n H[ 2 

H21 H22 


(116) 























where Hij = with i,j = 1, 2, ’s denote 

the basis states which are defined in Eq. (11141) . 

In Eq. d65j, there is a term (1 + h 2 )/2 on both left 
and right sides of the expression. We can let (1 + h 2 )/2 
operate on the wave functions of Hq and obtain new wave 
functions, then deal with the middle term of Eq. (l65l) with 
the new wave functions. 


1 

“ 2 


+ h 2 ) 


g(r) y™L 


m 2 ^ £ P 

UJ 2 ' Co! 2 

yp m2 I i 

Co! 2 Co> 2 ' 


\l(l 

\-2 V ( 


g(r) y™h 

y?,i B 

k ,—1 d 


- l )9( r ) + 2Z3j( V - Vj,L 
*[^(¥ + &9(r) + |(-^ + l)/(r)] vli 

//{'■) vTL 

; /b-) 


(117) 


the above equation shows the transformation from the 
wave functions of H 0 to the new wave functions with the 
operation of (1 + h 2 )/2. 

With the new wave functions the perturbative contri¬ 
butions of the H' to the eigenvalues of states can be cal¬ 
culated easily, they are given below: 

(1) J p = o- 

Only state with j = 1/2, l = 0 can construct the J p = 
0“ state, we have 


E n , o,|,o — En, o,| 

Vv ( 2 4o,| +<,0,i) 9n, 0,1 rdr 

(118) 

(2) J p = 1- 

Both states with j = 1/2, l = 0 and j = 3/2 ,1 = 2 can 
construct the J p = 1~ state, we have 



and 

E 2 3 1 = E 2 2 + 77 - 

n,2, 2 ,i n,2, 2 3mi 

-^,2,1/^2,f + 5r /n,2, 

(3) J P = 0+ 

Only state with j = 1/2, l = 
0+ state, we have 



C,2,f) rdr ( 12 °) 


1 can construct the J p = 


E n , i,±,o — En, i,| 

+ ^L Vv ( 29n>i> * +r9 ' n ’ i ’i) fn ^ rdr 

(121) 


(4) J p = 1+ 

Both states with j = 1/2, l = 1 and j = 3/2, Z = 1 can 
construct the J p = 1 + state, the matrix elements of the 
mass matrix are 




r 2 dr, 


( 122 ) 


+ 3 /n,i,f5„,i,i +r/ nili | 5 ; i i 

+^ n ,i 1 i/;,i i ! +r5 n ,i,|/; i i i i) rdr, 

#21 = #1*2, 



(123) 

(124) 


(125) 


With the matrix elements given above, one can get the 
eigenvalues of the two mixing states and the mixing angle 
easily by diagonalizing the mass matrix. 

(5) J p = 2+ 

Both states with j = 3/2, l = 1 and j = 5/2, l = 3 can 
construct the J p = 2 + state, we have 



7r/„,i, *g' n ,i ,§ 


(126) 


and 


^n, 3,§,2 ~^n, 3,-| 





(127) 


(6) J p = 2- 

Both states with j = 3/2, l = 2 and j = 5/2, l = 2 can 
construct the J p = 2~ state, the matrix elements of the 
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mass matrix are 


mass matrix are 


# n = E n 


H „ = E, 


1 r°° / 

+ 5wl’ J Q Vv “ 7r 3n,2,|/ rl ,2,| 

+ 3r /n, 2 ,|ffl, 2 ,|) rdr ! ( 128 ) 


#12 = 


5^V^l ^( 3 /n,2,|Sn,2,f 


+ 3 /n,2,§Jn,2,§ + r /ra,2,§S( lj 2,§ + r /n,2,§ #1,2,§ 

+rg n ,2,%fn 2 S + r 5n,2j/^ 2 a) rdr ’ ( 129 ) 

#21 = #1*2, (130) 

2 /, °° 


#22 = E„ 


V v 


( 42 /n, 


2jSn.2,# 


5mi Jo 

+ 4 ^n,2,|/^ 2 ,| - r /n,2,|<,2,|) r rfr - ( 131 ) 


By diagonalizing the mass matrix, one can get the eigen¬ 
values of the two mixing states and the mixing angle. 

(7) J p = 3“ 

Both states with j = 5/2, l = 2 and j = 7/2, l = 4 can 
construct the J p = 3“ state, we have 


n, 3 ,t 


+ 


7m 


1 Jo ^ r 3n,3,f /n,3,| 

+ 2 r ’/ra,3,|5n,3,|) r dr , (l 34 ) 

^ f V » ( 3 /»,3, 


#12 = 


7mi 


+ 3 /n,3,|5n,3,| +^/n,3,|Sn, 3 ,| + r /n,3,l5( l ,3,| 

+ r 5n,3,|/^, 3 ,} + r 5n,3,i/( 1 ,3 1 |) r dr, (135) 

(136) 

j /*oo 


#21 = # 1 * 2 , 

#22 = #„, 3 ,| - 


7m 


V v 


i Jo 


( 40 /n,3, 


|5n,3,| 


+ 11 r 9n,3,lfn, 3 l ~ 3 T fn, 3, f<4, 3) |) P dr - ( 137 ) 


The eigenvalues of the two mixing states and the mixing 
angle can be obtained by diagonalizing the mass matrix. 

Once the eigenfuncion of H 0 is obtained, the pertur¬ 
bative correction of #' can be calculated directly by fol¬ 
lowing the procedures presented above. 



3 r /n,2,#5 n .2.i 


(132) 


and 



llrfn,4,!<4,z. 

(133) 


(8) J p = 3+ 

Both states with j = 5/2, l = 3 and j = 7/2, Z = 3 can 
construct the J p = 3 + state, the matrix elements of the 


IV Numerical result and discussion 


The parameters used in this work are the quark masses, 
m u ,d ., m s , m c , mb and three potential parameters a s , 6, 
c. In the calculation of the spectra, we find that the 
model is capable to give reasonable spectral structure for 
each meson system. In the framework of Bethe-Salpeter 
equation, the parameter b is responsible for the energy 
levels of states with higher quantum numbers n or l. We 
find that if the parameter b is taken as a constant for 
all the meson systems, the energy levels for the excited 
meson states decrease as the value of m 2 decreases. As a 
result, if we determine the value of b for D s , B s mesons, 
then calculate D , B mesons with the same value of b 
but smaller value of m 2 , we can not get the correct en¬ 
ergy levels for the excited states of D, # mesons. This 
can be explained by the form of # 0 . Unlike the Dirac- 
like Hamiltonian, the influence of the confining potential 
V s (r) weakens as m 2 decreases in Eq. (l64l) . This diffi¬ 
culty suggests that the parameter b may depend on the 
masses of the quark and antiquark, especially the lighter 
one of them. 

With the considerations above, our fitting of the pa- 
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rameters gives the following values 

m u ,d = 0.360 GeV, 
m s = 0.550 GeV, 
m c = 1.478 GeV, 
m b = 4.865 GeV, 
a s = 0.513, 

f 0.350 GeV 2 for (cq, bq) system, 
\ 0.260 GeV 2 for (cs, bs) system, 
c = -0.320 GeV. 


Numerical calculation shows that the solution of the 
wave equation is stable when L > 5 fm, N > 35. Here 
we take L = 10 fm, N = 50. We fit the spectra of the 
heavy-light D, D Sl B , B s mesons mainly based on the 
meson states presented in PDG (24}. The numerical re¬ 
sults for the spectra of heavy-light mesons are presented 
in two tables. Table I is for D , D s mesons and Table II 
for B , B s mesons. The obtained spectra are in reasonable 
agreement with the experimental measurements. Theo¬ 
retical deviations from experimental data mainly appear 
in the D s meson sector, where the calculated mass of 
the _D* 0 (2317)° resonance is the worst. In the two ta¬ 
bles, one can see that apart from the D* 0 (2317)° and 
H s i(2460) resonances, mass calculations for the rest of 
the resonances are in good agreement with experimen¬ 
tal data. The discrepancy may be ascribed to the naive 
assumption of the kernel, the dropped a 2 term in the 
Hamiltonian can also make the result incomplete. The 
highly excited meson states are also calculated in the 
spectra and the newly observed charmed meson states 
are identified in our model. 

As for D mesons, several resonances are observed by 
LHCb collaboration in the mass region between 2500 and 
2800 MeV as well as thejregion around 3000 Mev, their 
masses are measured as @] 


M(Dj{ 2580)°) 
M{D*j{ 2650)°) 
M(Dj(2740)°) 
Af(D}(2760)°) 
M(D*j( 2760)+) 
M(Dj(3000)°) 
M(D}(3000)°) 
M(D*j{ 3000)+) 


2579.5 ±3.4 ±5.5 MeV 
2649.2 ± 3.5 ± 3.5 MeV 
2737.0 ±3.5 ±11.2 MeV 

2760.1 ± 1.1 ± 3.7 MeV 

2771.7 ±1.7 ±3.8 MeV 

2971.8 ± 8.7 MeV 

3008.1 ±4.0 MeV 

3008.1 MeV. 


The assignments of the above observed states are listed 
in Table I. In the calculated spectrum of D meson, 
Dj{ 2740) can be identified as the |1 5 / 2 D 2 ) state with 
J p = 2“ and 11}(2760) can be identified as the |1 5 / 2 D 3 ) 
state with J p = 3 _ , which agrees with Ref. E3 . In 
Ref. [§], Dj{ 2580) and _D}(2650) are identified as n = 2 
excited states with J p = 0“ and J p = 1”, respectively. 
Our result favors this assignment with Dj (2580) assigned 


as the |2 1 / 2 5'o) state and £>}(2650) as the \2 l t 2 S\) state. 
One can see that £>}(3000) and Dj( 3000) can be assigned 
in Table I as the |2 1,/2 Pi) and |2 3 / 2 P 2 ) states, respec¬ 
tively. But as presented in Ref. [8|, £>}(3000) favors the 
natural parity and Dj( 3000) the unnatural parity. With 
this consideration, D}(3000) is identified as the |2 3 / 2 P 2 ) 
state with J p = 2 + and Pj(3000) as the |1 7 / 2 P 3 ) state 
with J p = 3 + in this work. 

As for D s mesons, Several new D s meson states 
have been observed besides the low-lying states, such as 
P sJ (2632) 0 , (2860) @, D; 1 (2710) |, D sJ { 3040) 

4|. D s j{ 2632) is measured to be 2632.5±1.7 MeV in Ref. 
l|, we assign it as the |2 1 / 2 5'o) state with J p = 0 _ . For 
-0^(2710), it is measured to be 2708 ± 91^ Mev in Ref. 
0, it can be assigned as the |2 1 / 2 5 i) state with J p = 1~, 
which agrees with Refs. [lH UK]. The D s j(3040) reso¬ 
nance is observed by BABAR Collaboration at a mass 
of 3044 ± 8± 30 Mev 0. Several possible assignments of 
this resonance is discussed in Ref. [§]. It is identified as 
the |2 1 / 2 Pi) state with J p = 1+ in our predicted mass 
spectrum for D s meson. Recently, LHCb collaboration 
identifies the resonance D* sJ { 2860) as an admixture of 
spin-1 and spin-3 resonances and determines the masses 
and widths of the two states to be 


M(D* sl ( 2860)“) 

r(^(2860)-) 

M(D* s3 ( 2860)-) 
T(D* s3 ( 2860)-) 


2859 ± 12 ± 6 ± 23 MeV, 

159 ± 23 ± 27 ± 72 MeV, 
2860.5 ± 2.6 ± 2.5 ± 6.0 MeV, 
53 ± 7 ± 4 ± 6 MeV. 


In Table I, One can see that our predictions for the masses 
of the |1 3 / 2 Di) and \1 5 ^ 2 D 3 ) states are around 2860 Mev. 
Our model appears to be able to interpret the two states 
being the J p = 1~ and 3 _ members of the ID family. 
We can assign P* 1 (2860) _ as the |l 3 / 2 P>i) state with 
J p = 1~ and D* 3 (2860)“ as the |1 5 / 2 P 3 ) state with 
J p = 3". 

The wave function of each bound state can be obtained 
simultaneously when solving the wave equation. The ra¬ 
dial wave functions g n ,i,j(r) and f n ,l,j(r) for D meson are 
depicted in Fig.l as an example. [ g 2 ; .(r) ± / 2 ; j(r)]r 2 
is the possibility density distributed along the quark- 
antiquark distance r. In Ref. 37], the radial wave func¬ 
tions are also presented, they are calculated in a model 
derived by reducing the spectator equation in relation to 
the Bethe-Salpeter equation [§1]. One can find that their 
wave functions share the same pattern with ours, though 
the Hamiltonian Hq in our model looks very different 
from the Dirac-like Hamiltonian used by them. 


V Summary 

We construct a relativistic model for heavy-light 
quark-antiquark systems by studying the reduction of the 
instantaneous Bethe-Salpeter equation. The kernel we 
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TABLE I: Spectrum for P and D„ mesons. E° denotes the lowest order energies. P phys ' includes all the corrections of order 
1/rriQ. All units are in MeV. 


n J Lj 

Meson 

Pexpt. [8, J24| 

P u 

^phys. 

Meson 

Pexpt. [241 

P u 

j^ptiys. 

1 1/2 S 0 

D 

1869.62 ±0.15 

2105 

1859 

Df 

1968.49 ± 0.32 

2179 

1949 

l^Sl 

D * 

2010.28 ±0.13 

2105 

2026 

p:± 

2112.3 ±0.5 

2179 

2110 

1 1/2 P 0 

Dq (2400)° 

2318 ± 29 

2644 

2357 

P*o(2317)° 

2317.8 ±0.6 

2659 

2412 

l 1/2 Pi 



2644 

2529 

Psi (2536) 

2535.12 ±0.13 

2659 

2562 

l 3/2 Pi 

Pi (2420) 

2421.3 ±0.6 

2530 

2434 

P s i(2460) 

2459.6 ± 0.6 

2618 

2528 

1 3/2 P 2 

P 2 (2460) 

2464.4 ± 1.9 

2530 

2482 

P* 2 (2573) 

2571.9 ±0.8 

2618 

2575 

1 3 / 2 Pi 



2968 

2852 

P^ (2860) ^ 

2859 ± 12 ± 6 ± 23 [5] 

2976 

2873 

1 3/2 P>2 



2968 

2900 



2976 

2916 

1 5 ' 2 d 2 

Pj(2740)° 

2737.0 ±3.5 ±11.2 

2793 

2728 



2889 

2829 

1 5//2 P>3 

P}(2760)° 

2760.1 ± 1.1 ±3.7 

2793 

2753 

P 3 * 3 (2860)- 

2860.5 ± 2.6 ± 2.5 ± 6.0 [5] 

2889 

2852 

i 5/2 p 2 



3188 

3107 



3201 

3128 

1 5 / 2 P 3 



3188 

3134 



3201 

3152 

17/2^3 

Pj(3000)° 

2971.8 ±8.7 

2996 

2942 



3098 

3049 

2 1/2 S 0 

Pj(2580)° 

2579.5 ±3.4 ±5.5 

2821 

2575 

P s j(2632) 

2632.5 ± 1.7 [1] 

2849 

2624 

2 1 / 2 S 1 

P}(2650)° 

2649.2 ±3.5 ±3.5 

2821 

2686 

Pal (2710) 

2708 ± 9±[?, [31 

2849 

2729 

2 1 / 2 Po 



3133 

2902 



3136 

2918 

2 1/2 Pi 



3133 

2999 

P s j(3040) 

3044 ± 8+f [4] 

3136 

3017 

2 3/2 Pi 



3052 

2932 



3106 

2994 

2 3/2 P 2 

P}(3000)° 

3008.1 ±4.0 

3052 

2969 



3106 

3031 

2 3/2 P>i 



3350 

3228 



3360 

3247 

2 3/2 P> 2 



3350 

3260 



3360 

3278 

2 5/2 P> 2 



3225 

3139 



3296 

3217 

2 5/2 P 3 



3225 

3160 



3296 

3237 

2 5/2 P 2 



3516 

3425 



3534 

3449 

2 5/2 P 3 



3516 

3444 



3534 

3468 

2 7/2 P 3 



3371 

3301 



3455 

3390 


choose is based on scalar confining and vector Coulomb 
potentials, and it shows a Coulombic behavior at short 
distance and a linear confining behavior at long distance. 
The bound states of D , D S: B , B s mesons are studied in 
this model. The predictions of the spectra and the assign¬ 
ments for the newly observed charmed mesons are pre¬ 
sented. Our results are in reasonable agreement with ex¬ 
perimental measurements. The spectral structure is ob¬ 
tained in the relativistic model with merely three poten¬ 
tial parameters. We find that the parameter b in the con¬ 
finement potential depends on the masses of the quarks 
in the framework of Bethe-Salpeter equation. Theoreti¬ 
cal deviations from experimental data mainly appear in 
the D s meson sector, especially for the D* 0 (2317)° and 
D s i(2460) resonances. The discrepancy may be ascribed 
to the naive assumption of the kernel, the influence of 
other kernels with different spin structure can be studied 
in further researches. The wave functions of each bound 
state can be obtained by solving the wave equation and 
be used in the study of B and D decays. 
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Appendix A: The reduction for double-heavy 
systems 


It is interesting to discuss the differences of the reduc¬ 
tion results with and without the approximation which 
is mentioned in section II, 


Mp) -a 1 . 


(138) 


In this paper we discard the approximation and retain 
the original form of the Salpeter equation. Here we only 
consider systems with vector interaction for the sake of 
simplicity. 

As for the reduction of the double-heavy system, the 
procedure of the Foldy-Wouthuysen transformation has 
been generalized by Chraplyvy Hil in such a way that 
it can be applied to the two-body problem. If we write 
the general two-body Hamiltonian as 

H 12 = i + /3 (2 ) to 2 +££ + 80 + 08 + OO, (139) 
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TABLE II: Spectrum for B and B a mesons. E° denotes the lowest order energies. E phys ' includes all the corrections of order 
1/mQ. All units are in MeV. 


nP Lj 

Meson 

Eexpt. [241 

£ u 

^phys. 

Meson 

Eexpt. [241 


£ P hys. 

1 1/2 S 0 

B 

5279.25 ±0.17 

5362 

5262 

B s 

5366.77 ± 0.24 

5435 

5337 

l^Si 

B* 

5325.2 ± 0.4 

5362 

5330 

b: 

5415.4±| 

5435 

5405 

l 1/2 Po 



5859 

5740 



5883 

5776 

l 1/2 Pi 



5859 

5812 



5883 

5841 

l 3/2 Pi 

Bi(5721) 

5723.5 ± 2.0 

5772 

5736 

B s i(5830) 

5829.4 ± 0.7 

5860 

5824 

1 3/2 P 2 

B 2 (5747) 

5743 ± 5 

5772 

5754 

#2 (5840) 

5839.7 ± 0.6 

5860 

5843 

1 3//2 Di 



6174 

6128 



6188 

6146 

i 3/2 d 2 



6174 

6147 



6188 

6163 

i 5/2 d 2 



6013 

5989 



6108 

6085 

i 5/2 d 3 



6013 

5998 



6108 

6094 

i 5/2 f 2 



6375 

6344 



6392 

6363 

1 5/2 e 3 



6375 

6354 



6392 

6373 

i 7/2 f 3 



6195 

6175 



6295 

6276 

2 1/2 S 0 



6013 

5915 



6052 

5961 

2 1 '' 2 Si 



6013 

5959 



6052 

6003 

2 1/2 P 0 



6302 

6211 



6315 

6227 

2 1/2 Pi 



6302 

6249 



6315 

6266 

2 3/2 Pi 



6232 

6189 



6293 

6249 

2 3/2 P 2 



6232 

6200 



6293 

6263 

2 3/2 Pi 



6506 

6458 



6523 

6478 

2 3/2 P 2 



6506 

6471 



6523 

6491 

2 5/2 D 2 



6389 

6357 



6465 

6434 

2 5/2 P 3 



6389 

6365 



6465 

6441 

2 5/2 P 2 



6656 

6621 



6680 

6647 

2 5/2 P 3 



6656 

6629 



6680 

6654 

2 7/2 P 3 



6519 

6493 



6605 

6580 


then the transformed Hamiltonian has the form 


sign can be rewritten as: 


H 12 = Uy 2 l HUy 2 

= /3 (1) mi + p (2) m 2 +££+ ^(0£) 2 + ^-{£Of 

lull 2777-2 

/3 (1) /3 (2) 1 


+ 


8m? 

1 

8 TO 2 


(0£T - (Mr + Tr-s[[0£,££],oe\ 


8rri2 

[[£ 0 ,££\,£ 0 \ 


8 m{ 


pWp™ 

4miW2 


{{ 0 £, 00 },£ 0 } 


pW + p& 

4 (mi + m 2 ) 


(OO ) 2 


(140) 


E / E / ■ p\ 

2 \ u>i ) 2 \ wj j 

+awi +aw 2 + — {hi + h 2 )U— (h\ + h 2 ) 

= -y£ (1) ~ jP (2) +££ + 0£+£0 + 00, (143) 


Eq. (1151) can be written as (the wave function in the where U is defined as 
equation is omitted): 


-~^E - yE +AW 1 +AW 2 + 2 (fri + h 2 )U-(hi + h 2 ) 

= —mi - m 2 , (141) 

where 


U = Vi + V 2 , (144) 

Vi = V, (145) 

V 2 = — — [cc^ • a/ 2 ) + (a* 1 ' • r)(a^ ■ r)]Vi(146) 


A w» = a;* - mf, * = 1,2 (142) 

now we consider performing the Foldy-Wouthuysen 
transformation on Eq. (11411) . the left side of the equal 


From Ea. (1143jl . we can obtain ££, 0 £, £0 and OO. 
Inserting the four operators into Eq. (11401) and expanding 
the equation to order (1/E) 2 , we have the transformed 
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FIG. 1: The radial wave functions gn,i,j(r ) and for D meson as an example. The wave functions are defined in Eq. 

(TTOTl) . They are the radial part of the solution of the eigenequation associated with Hq. 


Hamiltonian 


Hi2 = UJl + w 2 + V — 


1 


4?71iTO2 


{a 


( 2 ) . 


^[[«" p.H« l2) -p] + i[l 

/3W+/3< 2 ) 


8 TO 2 


a 


4Eq 


-{ooy. 


p,{a {1) p,V 2 }} 
(1) p, V],a (1) -p\ 

(147) 


All the terms of the transformed Hamiltonian are the 
same as that of the Breit interaction except for the last 
term. In our scheme, we have 


00= ItfW +/3W)V 2 ^(I3W + 0W), 


/gw + pw 

4£’o 


{(DO) 2 -l 


Yt 

2 Eq' 


(148) 

(149) 


in the heavy quark limit, Eq ss m\ +m 2 - While with the 
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approximation, the result is 


OO = V 2 , 


g (1) +^ ) (oof ^ _vi _, 

4(?m+m 2 ) 2(mi+m2)’ 


(150) 

(151) 


Eg. (11491) differs from Eg. (11511) by a minus sign. In con¬ 
clusion, for double-heavy system we can take the approx¬ 
imation and have a simpler form of wave equation if we 
disregard the last term in Eq. (11471) . 
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